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ABSTRACT We calculated the incidence, mortality, and case fatality rates for Caucasians and non-Caucasians during 19th century 
yellow fever (YF) epidemics in the United States and determined statistical significance for differences in the rates in different 
populations. We evaluated nongenetic host factors, including socioeconomic, environmental, cultural, demographic, and ac- 
quired immunity status that could have influenced these differences. While differences in incidence rates were not significant 
between Caucasians and non-Caucasians, differences in mortality and case fatality rates were statistically significant for all epi- 
demics tested (P < 0.01). Caucasians diagnosed with YF were 6.8 times more likely to succumb than non-Caucasians with the 
disease. No other major causes of death during the 19th century demonstrated a similar mortality skew toward Caucasians. Non- 
genetic host factors were examined and could not explain these large differences. We propose that the remarkably lower case 
mortaUty rates for individuals of non-Caucasian ancestry is the result of human genetic variation in loci encoding innate im- 
mune mediators. 

IMPORTANCE Different degrees of severity of yellow fever have been observed across diverse populations, but this study is the 
first to demonstrate a statistically significant association between ancestry and the outcome of yellow fever (YF). With the global 
burden of mosquito-borne flaviviral infections, such as YF and dengue, on the rise, identifying and characterizing host factors 
could prove pivotal in the prevention of epidemics and the development of effective treatments. 
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The identification of host factors that aid or restrict pathogens 
has transformed our understanding of disease susceptibility, 
outcomes, and possible treatment responses for major global pub- 
lic health threats, including malaria (1) and HIV/AIDS (2). Un- 
derstanding yellow fever (YF) host factors is a crucial part of a 
comprehensive approach to protect the public, as information 
about host variation, when combined with climate, geographical, 
and additional population data, could be used to better predict the 
future locations of especially severe outbreaks of YF. Most impor- 
tantly, an understanding of critical host factors could help develop 
novel antiflaviviral treatments. 

Identification of host factors for mosquito-transmitted flavivi- 
ruses, important global and reemerging threats, has begun using 
cell-based assays (3-5), but studies of human populations are only 
beginning. In fact, genome-wide association studies (GWAS) have 
recently been used to identify potential genetic host factors for 
dengue hemorrhagic fever (DHF)/dengue shock syndrome (DSS) 
(6) and also nonsevere dengue fever (DF) (7). Estimates typically 
put the global burden of YF at 200,000 cases and 30,000 deaths 
annually (8), although a recent study suggests that these numbers 
may significantly underestimate this burden (4). The reemergence 



of YF has been attributed to lack of vaccination in regions at risk 
for YF, the spread of the mosquito vector range, and increased 
urbanization (9). These factors, in tandem with increased global- 
ization and climate change, contribute to a growing number of 
areas at risk for YF; of particular concern is the potential of an 
urban epidemic in a tropical or subtropical megacity (9). 

The United States in the 19th century was a unique setting to 
observe differences in YF susceptibility between populations of 
different ancestry, as the country contained large populations of 
European and African ancestries living in close proximity. Fur- 
thermore, the United States was never an area where YF was en- 
demic, and therefore, few people contracted YF outside of re- 
corded epidemics, decreasing the ambiguity of diagnosis, which 
complicates historical evaluations of disease variance between 
populations. Anecdotal evidence from 19th century YF epidemics 
in the United States suggests that Caucasians succumbed to the 
disease at considerably higher rates than those of African descent 
(10), but this difference in mortality has never been rigorously 
established. Studies on historical epidemics (10, 11) and research 
focusing on responses to the YF vaccine (12) have used data trends 
to suggest that differences in mortality between races could be 
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TABLE 1 Case fatality rates stratified by race during yellow fever epidemics in the United States (1808 to 1878)" 
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explained by genetic factors. This explanation and other explana- 
tions for the different susceptibilities (13) have not been subjected 
to statistics-based hypothesis testing. 

In this study, we integrated data from a variety of sources, 
including primary-source documents created by medical person- 
nel who treated patients during the epidemics, census data, and 
information about the epidemics compiled by contemporary his- 
torians. Taking advantage of unique data sets from 19th century 
U.S. urban YF epidemics, we utilized tests of statistical inference to 
establish significant differences in YF mortality and case fatality 
rates between Caucasians and non-Caucasians. By performing 
rigorous hypothesis testing on the effects of socioeconomic, envi- 
ronmental, cultural, and demographic factors on this association, 
our results strongly suggest that host genetic variation affects yel- 
low fever virus (YFV) infection outcomes. Based on these results, 
a large-scale genetics study to identifj? host factors for YF suscep- 
tibility is warranted. 

RESULTS AND DISCUSSION 

Caucasians were more likely to succumb to yellow fever than 
non-Caucasians during 19th century epidemics. There was no 
consistent difference in the incidence of YF between Caucasians 
and non-Caucasians (see Text SI and Table SI in the supplemen- 
tal material). On the other hand, mortality rates, which varied 
widely between epidemics, were 1.86- to 39.42-fold higher for 
Caucasians in all 16 epidemics for which data are available (Ta- 
ble S2). The differences in mortality rates between populations 
were statistically significant {P < 0.01), indicating that 
population-based differences in YF susceptibility existed. 

Since the set of clinical features of severe YF cases is distinctive, 
the potential of mistaking the cause of death as YF was unlikely. 
Nonetheless, it is possible that mortality rates could have been 
significantly affected by differences in socioeconomic status, par- 
ticularly access to health care providers. Cultural bias could also 
have affected mortality rates, since until the 1840s many physi- 



cians believed that those of non-Caucasian descent could not con- 
tract YF (14), a misconception which has been refuted. To mini- 
mize these confounding variables, we focused on epidemics where 
we could analyze case fatality rates, which utilize only the cases 
diagnosed by medical personnel. Six epidemics met this additional 
inclusion criterion. 

In these six epidemics, the case fatality rate for Caucasians var- 
ied from 25.0 to 72.5%, while the case fatality rate for non- 
Caucasians varied from 1.1 to 14.1% (Table 1). This remarkable 
skew in case fatality rates between races was observed in all six 
epidemics. When we constructed regression models for each of 
the six epidemics, the increased case fatality rates for Caucasians 
relative to non-Caucasians were all significant (P < 0.01 [Ta- 
ble 1]). When the data sets were combined, the case fatality rate 
was 6.8 times greater for Caucasians than non-Caucasians, and the 
odds of Caucasian patients succumbing to YF were 14.6 times 
greater than those of non-Caucasian patients (with a 95% confi- 
dence interval of 13.59 to 15.76 [Table 1]). 

We considered recording bias as an explanation for this dra- 
matic difference in case fatality rates. For instance, it is possible 
that non-Caucasians were poorly followed up by medical person- 
nel. This confounding variable was minimized for the epidemic 
among the U.S. troops in New Orleans in 1867, given the War 
Department's emphasis on record-keeping (15). The Caucasian 
troop case fatality rate was 1.99-fold higher {P < 0.01) and the 
odds of Caucasian troops dying from YF were 2.56 times greater 
than those of non-Caucasian soldiers (Table 1). Given the more- 
homogeneous treatment among troops, these data also argue 
against the possibility that YF "treatments" such as quinine, castor 
oil, and calomel (15, 16) or some other aspect of the care (17) 
could have increased case fatality among Caucasians, a population 
that was more likely to have access to care, particularly earlier in 
the disease cycle or for milder cases (see Text S2 in the supplemen- 
tal material). We conclude from these data that there was a re- 
markably higher and statistically significant case fatality rate 
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FIG 1 Time course of the YF epidemic in Memphis, TN, in 1878 was similar 
in Caucasians and non-Caucasians. The numbers of deaths from yellow fever 
over time during the epidemic in Memphis, TN, in 1878 are shown. 



among Caucasians compared to individuals of African descent 
during 19th century YF epidemics in the United States. 

Testing potential explanations for population-based differ- 
ences in YF mortality. Biological, demographic, cultural, socio- 
economic, or environmental factors may have contributed to the 
population-based differences in YF severity, so existing data were 
analyzed to assess whether each of these factors could have signif- 
icantly affected these rates. 

(i) It is unlikely that differences in viral genotypes explain 
differential severity within epidemics. As can be noted from Ta- 
ble 1, and consistent with what is known about dengue (5), sever- 
ity can vary between epidemics, and this could be attributed in 
part to the differential virulence of viral genotypes. Since the ad- 
dresses and workplaces of Caucasians and non-Caucasians within 
a city were often segregated ( 14), we considered the possibility that 
the differences in severity within epidemics were because the two 
races were being affected by two different strains of YF V. Given the 
unavailability of YFV isolates from these epidemics, we examined 
historical data to ascertain whether there was evidence for differ- 
ences between outbreaks among Caucasians and non-Caucasians. 
We investigated the kinetics of the YF epidemic in Memphis, TN, 
in the two populations by charting deaths from YF over time and 
found similar kinetics (Fig. 1). The consistency in disease kinetics 
between populations was mirrored when examining cases diag- 
nosed and deaths over time in the 1867 YF epidemic in New Or- 
leans, LA, among U.S. Army troops (see Fig. SIA and SIB in the 
supplemental material). The time courses for YF-associated 
deaths were similar in both populations, which is most consistent 
with one YF genotype spreading in both populations. Moreover, 
both the Memphis 1878 and the U.S. Army 1867 data indicate that 
at all times during an epidemic, the case fatality rate for Cauca- 
sians greatly exceeded that for non-Caucasians. 

(ii) Acquired immunity is unlikely to explain differences in 
yellow fever case fataUty rates. Population-based differences in 
YF outcome have been attributed to differing rates of acquired 
immunity. Immigration from Africa early in the 19th century, 



immigration of Caucasians from regions where YF is not endemic 
later in the century, and different behavioral or cultural patterns 
between the two groups could contribute to differing levels of ac- 
quired immunity between the two populations (13). It is important 
to note, however, that individuals who had previously acquired 
immunity were not included in the calculation of case fatality 
rates, because multiple clinical infections of YF in one individual 
have not been documented in the historical literature (15, 18) 
and are considered exceedingly rare (http://www.cdc.gov/travel 
-training/local/HistoryEpidemiologyandVaccination/page24492 
.html). Therefore, for this important reason, and others discussed 
in detail in Text S3 in the supplemental material, we conclude that 
acquired immunity was not a determinant of differential YF case 
fatality between Caucasians and non-Caucasians. The confound- 
ing issue of cross-reacting partial immunity due to prior infection 
with dengue viruses cannot be completely ruled out. Nevertheless, 
we argue against a major role for cross-reacting immunity for two 
reasons. First, it is likely that incidence of dengue fever, which was 
also epidemic in the 19th century United States and was transmit- 
ted by the same mosquito vector, was similar for Caucasians and 
non-Caucasians, as we note for YF (Table SI). Therefore, immu- 
nity levels should have been comparable for the two populations. 
Second, given what is known about cross protection between dif- 
ferent dengue serotypes, it is unlikely that these viruses could pro- 
vide lasting protection against YFV. 

(iii) Demographic factors are unlikely to explain differences 
in case fataUty rates for YF. Both age and gender could affect 
susceptibility to YF, indeed males may be more likely to succumb 
to YF than females (12) and children may have milder YF cases 
than adults (17). We examined whether demographic differences 
between Caucasian and non-Caucasian populations could explain 
the difference in YF case fatality. Case fatality rates stratified by age 
and population were analyzed (Table 2; see Table S3 in the sup- 
plemental material). Even when controlling for age and gender 
(there were only males in the data set), the YF case fatality rates 
were at least 1.86-fold higher for Caucasians (P < 0.01 [Table 2] ). 
In contrast, these population-based trends were not seen with 
cholera, as the differences in case fatality rates were the same or 
minimally lower in Caucasian than non-Caucasian soldiers (Ta- 
ble 2). 

To examine whether distortions in the underlying population 
structure could have affected skews within or between popula- 
tions, a case study on the 1878 YF epidemic in Memphis, TN, was 
performed. Although there were 13% more Caucasian males than 
females versus 18% more non-Caucasian females than males re- 
siding in Memphis in 1870 (see Fig. S2A and S2B and Table S4 in 
the supplemental material), these small distortions are not suffi- 
cient to explain a 7.44-fold difference in case fatality rate for the 
1878 Memphis YF epidemic (Table 1). Indeed, the U.S. Army data 
presented above (Table 2 and Table S3) demonstrated that among 
age-matched men, YF case fatality was greater for Caucasians. 
Therefore, while age and gender may have played a small role in 
the differences in YF case fatality rate between the two popula- 
tions, these factors could not account for the remarkable differ- 
ences observed. 

(iv) Socioeconomic, environmental, and cultural factors are 
unlikely to account for differences observed. Socioeconomic sta- 
tus could have had widespread ramifications for disease response, 
as it could have affected overall health, access to health care, and 
diet. The absence of any evidence for nutritional deprivation lead- 
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TABLE 2 Case fatality rates from yellow fever and cholera among U.S. Army troops in 1867' 



Disease and age 
of soldier'' 


No. of deaths/100 
cases of disease 
in Caucasians 
(total no. of cases) 


No. of deaths/100 
cases of disease 
in non-Caucasians 
(total no. of cases) 


Fold change in 
case fatality 
rates (C/nC)'^ 


Odds (95% confidence 
interval) of Caucasians 
versus non-Caucasians 
succumbing to disease 


Yellow fever 










Less than 20 yr 


30.0 (237) 


12.1 (33) 


5.67*** 


15.97 (6.02-55.30) 


20-29 yr 


30.2 (770) 


16.3 (123) 


1.86** 


2.23 (1.38-3.79) 


Cholera 










Less than 20 yr 


22.2 (45) 


50.0 (6) 


0.44 


0.29 (0.05-1.75) 


20-29 yr 


37.8 (185) 


39.1 (87) 


0.97 


0.95 (0.56-1.61) 



" Data from reference 15. 

^ The average age of soldiers less than 20 years old was not given in the primary source. 

The case fatality rate in Caucasians (C) to the case fatality rate in non-Caucasians (nC) is shown. Race was found to be a significant predictor of outcome in the logistic regression 
model as follows: **, P < 0.01; P < 0.001. 



ing to better outcomes and Civil War data for U.S. Army soldiers, 
which partially homogenized socioeconomic differences, suggest 
that lower socioeconomic status would not be expected to lead to 
improved YF outcomes. Eighty-one causes of death among U.S. 
Army soldiers in 1 864 were analyzed, and indeed, contrary to what 
was observed for YF, most causes of death were skewed toward 
non-Caucasians (Fig. 2; see Table S5 in the supplemental mate- 
rial). In fact, most other infectious causes of death exhibited a 
higher case fatality among non-Caucasians (e.g., typhus fever) 
(Fig. 2; see also discussion in Text S2 in the supplemental mate- 
rial). We note that YF is more distorted than all other infectious 
diseases and all causes of death among U.S. Army soldiers, except 
a category labeled "Other accidents and injuries," which likely 
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FIG 2 The yellow fever case fatality rates among U.S. Union troops in 1864 
demonstrate a skew toward Caucasians. The ratios of case fatality rates for 
Caucasians over the case fatality rate for non-Caucasians for 8 1 selected causes 
of death among U.S. Union troops in 1864 are shown. The total numhers of 
deaths from different causes follow: 223 from typhus fever, 64 from mumps, 
1,445 from remittent fever, 1,922 from measles, 1,917 from typho-malarial 
fever, 223 from typhoid fever, 3,031 from smallpox and varioloid, and 336 
from yellow fever. 



reflects the fact that non-Caucasians were deployed in battle less 
frequently than Caucasians (Fig. 2). We conclude from these com- 
parisons that overall health or health care disparities are very un- 
likely to explain the skew toward higher case fatality rate for Cau- 
casians. Furthermore, our analysis indicates that cultural factors 
did not highly impact this skew (Text S2). 

The weakening of alternative hypotheses, which invoke non- 
genetic factors, led us to conclude that genetic factors explain the 
different YF responses in Caucasians and non-Caucasians. As a 
result, genetic differences, particularly those contributing to in- 
trinsic and innate immunity, should be directly explored. This is 
fully consistent with proposals that populations whose ancestors 
lived in areas where YF is endemic would be more likely to have 
undergone selection for resistance that increase YF survival (10), 
similar to the malarial resistance afforded to those of African de- 
scent (1). While the research presented cannot formally confirm 
the possibility of host genetic differences significantly influencing 
YF severity, it is an important first step toward this goal. It should 
be noted that before allelic variation in IL28B was implicated in 
treatment-induced hepatitis C virus clearance, genetic variation 
was inferred based on strong differences between individuals of 
dissimilar ancestry (19). By analogy, we posit that our work 
strongly suggests the need for future experiments that compare 
genotypes of individuals with different YF outcomes. 

MATERIALS AND METHODS 

Study design and definitions. Reviews listing 19th century YF epidemics 
(11,20,21) were used to identify periods labeled at the time as "YF epi- 
demics" by medical personnel or organizations (20) and to locate 
primary-source documents about these epidemics. The case numbers or 
numbers of deaths from YF by population were sought in reading these 
primary sources. If these data were available, then the epidemic was sub- 
jected to rigorous selection criteria, which are defined immediately below, 
to determine whether it would be included in this research. 

To be included, an epidemic must have occurred during the 19th 
century and labeled as such by medical personnel at the time. We re- 
stricted the analysis to epidemics consisting of at least 50 documented 
cases of YF within a 5-month period (midsummer to fall of the same year) 
in the same geographical area. The case number cutoff was determined to 
mitigate the effects of misdiagnosis, and the time period reflected when 
the climate would most likely support the mosquitoes, which transmitted 
YF. When clinical descriptions from an epidemic were available, the signs 
and symptoms were compared to the clinical features of YF recognized 
today (22). While some features of YF are nonspecific, such as fever, head- 
ache, myalgia and weakness, lumbosacral pain, anorexia, nausea, vomit- 
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ing, and a decreased pulse, the clinical signs of severe cases, jaundice, 
hematemesis (sometimes referred to as "black vomit" in the historical 
record), oliguria, bleeding from orifices, and delirium, when observed 
together and in the context of an epidemic, are considered pathogno- 
monic for YF (22-24). Sixteen out of 333 19th century U.S. YF epidemics 
(20) met these criteria. 

Study population. In each epidemic, an individual's race was used as a 
proxy for ancestry, Caucasian or non-Caucasian. Indeed, U.S. Census 
data reveal that over 97% of non-Caucasians living in the United States 
during the late 19th century were of individuals of African descent (25) 
(see supplemental material). 

Statistical analysis. Data analysis was performed in R Studio using the 
contingency tables for the raw data (http://www.rstudio.com/). Fold dif- 
ferences between Caucasians and non-Caucasians were found by dividing 
the Caucasian case fatality rate by the non-Caucasian case fatality rate. 

For each epidemic, a logistic regression of race on outcome was found 
in R Studio using Caucasians as the reference level. The odds ratios and the 
corresponding 95% confidence intervals were then found from the model. 
All the data sets for which an incidence rate, case mortality rate, or case 
fatality rate could be found were then combined, and three regression 
models were fit for the three types of rates. These models were analyzed in 
the same way as the individual epidemics. 

SUPPLEMENTAL MATERIAL 

Supplemental material for this article may be found at http://mbio.asm.org/ 
lookup/suppl/doi:10.1128/mBio.01253-14/-/DCSupplemental. 

Text 81, DOCX file, 0.1MB. 

TextS2, DOCX file, 0.1MB. 

Text S3, DOCX file, 0.1MB. 

Text S4, DOCX file, 0.2 MB. 

TextSS, DOCX file, 0.1MB. 

Figure SI, DOCX file, 1.2 MB. 

Figure S2, DOCX file, 1.5 MB. 

Table SI, DOCX file, 0.1 MB. 

Table S2, DOCX file, 0.1 MB. 

Table S3, DOCX file, 0.1 MB. 

Table S4, DOCX file, 0.1 MB. 

Table S5, DOCX file, 0.1 MB. 

Table S6, DOCX file, 0.1 MB. 
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